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1. Introduction 

In recent years proposals for quantum computation and spintronics have initiated search 
for semiconductor nonmagnetic bulk materials and heterostructures that might be used 
to realize high values of electron spin polarization, on the one hand, and its long storage 
time, on the other. At low temperatures, this task is successfully realized by exploiting 
low spin-relaxation rates of charge carriers localized in semiconductor quantum wells or 
confined in self-organized quantum dots [U EJ EJ H]. With increasing temperature up 
to room temperature the free carriers get delocalized and, in zinc-blende-lattice bulk 
semiconductors, their spin relaxation time drastically decreases to values of the order of 
100 ps or even less [31 El EE [7] which limits the electron spin polarization under continuous- 
wave (cw) photoexcitation to few percent [31 [5] • The situation can be partially improved 
in strongly n-doped samples where the intense electron-electron collisions essentially 
suppress the D'yakonov-Perel' mechanism of spin relaxation, as shown theoretically [8] 
and experimentally [9], see also a comprehensive review article by Wu et al [10] 

Recent experimental studies of electron spin dynamics in GaAsN alloys have 
revealed strong contradiction with the existing expectations: extremely high spin 
polarization, up to 90%, of the free electrons and its preservation during ~lns have 
been found at room temperature in nonmagnetic undoped GaAsN bulk films and 
GalnAsN quantum wells [UJ H21 [131 [TH [151 UHl UHl HH] - In these experiments 
electron spin polarization was created by optical pumping and detected by polarized 
photoluminescence (PL). It has been established that the anomalously high polarization 
and spin memory arise due to spin- dependent recombination of free electrons with 
deep paramagnetic centers present in the nitrogen alloys and subsequent dynamic 
spin polarization of the centers [131 El EE EE]. The latter acts as a spin filter 
preventing the spin majority photoelectrons to recombine and promoting the spin 
minority photoelectrons to disappear from the conduction band. As a result the total 
density of photoelectrons increases superlinearly with the increasing pumping power 
which leads to an enhancement of the edge PL intensity (up to a factor of 8) [13j EUJ El] 
as well as of room-temperature photoconductivity [22j [23]. Recent measurements of 
optically detected magnetic resonance (ODMR) [21j [21] and optical orientation of 
electron spins [131 [T8] have provided a proof that the paramagnetic centers in question 
are Gai self-interstitial defects and established the critical dependence of the defect 
concentration on the N composition, Ga(In)AsN alloy growth conditions and post- 
growth treatment, e.g., thermal annealing. 

The spin-dependent recombination (SDR) gives rise to formation of a coupled spin 
system of free and bound electrons where the free-electron polarization is intricately 
related to the spin state of paramagnetic centers. For the quantitative description 
of this coupled nonlinear system we proposed the SDR model p31 EE] based on the 
model introduced by Weisbuch and Lampel [25] and used later by Paget [26] in a 
qualitative analysis of low-temperature optical orientation of electron spins in GaAlAs 
and GaAs, respectively. Our model contains nine parameters. At present most of them 
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have been found for GaAsN experimentally. Particularly, the measurements of cw and 
time-resolved polarized luminescence under circularly-polarized photoexcitation in the 
presence of a transverse magnetic field have allowed to estimate the spin-relaxation 
times of free and bound electrons, lifetime of photoholes, absolute value of free-electron 
g factor and signs of the free- and bound-electron ^-factors. The g factor of bound 
electrons is determined from ODMR studies [23]. 

In this paper we present a systematic theoretical study of the SDR in the three- 
and two-charge-state models with a focus concentrated on the conduction photoelectron 
spin polarization and circular polarization of edge PL induced by circularly-polarized 
radiation. It is shown that the two-charge-state theory is compatible with available 
experimental data on GaAsN alloys. The influence and specific role of each model 
parameters is analyzed and a complete set of them is determined for a particular GaAsN 
alloy sample. The experimental curves are compared not only with the results of exact 
model calculation but also with approximate analytical equations derived in the limit 
of weak or strong pumping powers and small or high magnetic fields. 

The paper is organized as follows. In section [2J we extend the theory of 
Shockley-Read recombination to allow for the spin polarization and spin-dependent 
recombination of photocarriers. We start from a more general three-charge-state model 
of the SDR and then reduce the consideration to the two-charge-state model for the 
particular case of cw photoexcitation. The exact equations are derived for the PL 
intensity and circular polarization, and analytical equations are derived in the simple 
limiting cases. In section [31 the values of nine parameters of the two-charge-state 
model are estimated by using available experimental data and reasonable theoretical 
considerations. Experimental and theoretical dependencies of the SDR ratio and degree 
of PL circular polarization on the pumping intensity and magnetic field are presented, 
compared and discussed in section [31 Section [5] contains the concluding remarks. 

2. Theory 

2.1. Three- charge- state model 

In solids, impurity atoms can form several energy levels in the band gap. To illustrate 
we refer to four-level Au impurity in germanium [27] and three-charged states of 
interstitial boron in silicon [28] or Fe, Rh in photorefractive crystals [291 E3- The 
multicharged states of Ga interstitials in GaAs are analyzed in reference [31], see 
also references therein. The statistics of defects with several trapping levels in 
semiconductors was analyzed decades ago [32] . Until now the rate equations describing 
light-induced changes in electron and hole densities and occupation of three-charged 
impurity states has been treated for spin-unpolarized charge carriers [291 EO]- In this 
section we generalize the consideration to the spin polarized photoelectrons and their 
spin-dependent recombination via deep paramagnetic centers. 

The three-charge-state model is summarized in figure 1(a). We assume that each 
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Figure 1. Band diagrams of the three-charge-state (a) and two-charge-state (b) 
models. Upper arrows indicate the excitation (up arrows) and recombination (down 
arrows) of conduction electrons described by the coefficients e„,c„ for the (2/1) 
transition between the two- and one-electron defect states and e' n ,c' n for the (1/0) 
transition. Lower arrows illustrate the hole excitation (e p ,e' p ) and recombination 
(c p , c p ). 



deep center can be in one of three states Di (i = 0,1,2) representing the center 
containing, respectively, no electrons, one bound electron with uncompensated spin 
±1/2 and two bound electrons in the singlet state with the zero total spin. The 
corresponding densities are denoted as N ,Ni,N2 and satisfy the condition 

N + N 1 + N 2 = N C , (1) 

where iV c is the total density of the deep-level centers. 

Under normal incidence of the circularly polarized light in the absence of a magnetic 
field, the electronic spins are polarized along the excitation direction (the axis z). Let us 
introduce the notations n±, N±, respectively, for the concentration of the free electrons 
and paramagnetic centers with the electron-spin component ±1/2 of the electron spin. 
Note that N + ± iV_ = N\. Hereafter, we assume the spin relaxation of photoholes to 
be very fast, neglect their spin polarization and use the notation p for the hole density. 
Taking into account interband optical generation of free carriers, their recombination 
and thermal ionization from impurity levels, the rate equations for the electron densities 
n±, the hole density p and the densities N±, N at zero magnetic field are 

dn± + n ± -n ± + ^ = + _ + + ^ ^ 

at It s 

dN± N± - N T 
+ 



dt 2r ; 



sc 

— -Rn.T ± Gn.± + R„ 4- — G„ _l + R v .+ — (?„.+ — R„ j- ± G' 



d ± + in±P = -R p -R p + G p + G' p + G opt , 
— -R'n + G' n + R' p — G' p . 
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Here are the photo-generation rates of free electrons with the spin ±1/2, G opt is 
the total generation rate G+ pt + G° pt proportional to the incident light intensity J, 7 is 
the coefficient for band-to-band recombination, the times r s and t sc describe the spin 
relaxation of free and single bound electrons (apparently, at room temperature r s <C r sc ). 
The rates for electron or hole recombination and thermal excitation are given by 

R n>± = 2c n N T n± , i4,± = c 'n N on± , 
R' n = R' n& + R' n> _ = c' n N (n + + n_) , 

R p = 2R Pj ± = c p N 2 p , R' p >± = c' p N±p , R' p = c' p (N + + N^)p , 

/ / / / / 

^n,± = ^"^2 , G n ± = e n iV ± , G n = G n+ + G n _ , 

G P)± = e p N ± , G p = e p (N + + AL) , G; = 2& Pt± = e' p N , 

the coefficients c n , e n , c p , e p describe the (2/1) processes Dx+e~ -H- D 2 and D 2 + h + <H- Z^x 
and have the same meaning as the corresponding coefficients in the original article by 
Shockley and Read [33] (see also the book [M]), the similar coefficients c' n , e' n , c' p , e' p refer 
to the (1/0) processes Do + e~ -H- D\ and Di + h + -h- where e~ and /i + symbolize 
an electron and a hole. 

In addition to equation (TjQ) the densities in consideration satisfy the neutrality 
equation 

p = n + N 2 - N - (iV 2 ° - iV °) , (3) 

where N® are the equilibrium densities given by [32] 

iV 2 ° 1 f n-E(2/l) \ N° (ii-E(l/0) 
exp — — , — = 2 exp 1 



JV? 2 fcsT y iV ° *\ k B T 

and Nq+Ni+N® = N c , where /j, is the electron equilibrium chemical potential, E(i+l/i) 
is the change in energy associated with adding one electron to the center in the process 
Di + e~ — > Di + i and the factors 2 and 1/2 take into account the double spin degeneracy 
of the paramagnetic state. 

The set of equations (DO)-© contains a lot of parameters. To simplify the analysis 
we assume that the photoexcitation rate G opt is high enough to neglect the thermal 
excitation rates G n> ±, G' ±, G p and G' . Moreover, we neglect the influence of band-to- 
band recombination on the kinetics of charge carriers, set 7 to zero in equations (j2J) 
and use the thus obtained values of n + ,n_ and p to calculate the photoluminescence 
(PL) intensity and polarization. For the sake of distinctness, we also assume that, in 
equilibrium, all the centers bind one electron per center, i.e., are paramagnetic and 
iV° = N c . 

Instead of four variables n±,N± we introduce the total free-electron density 
n = n + + n_ and spin S z = {n + — n_)/2, the density of paramagnetic centers (or 
centers with uncompensated spins) N\ = N + + N_ and the bound-electron total spin 
S C z = (N + — AL)/2. In a transverse magnetic field the electron spin exhibits the Larmor 
precession and, instead of the scalars S z and S cz , one should define the pseudovectors S 
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and S c pointing in the directions of the average free- and bound-electron spins. Then 
the set of kinetic equations (j2J) in the presence of magnetic field are rewritten as 
dn 

— + c n (N in - 4SS C ) + c' n N n = G opt , (4a) 

^ + {c p N 2 + c'pN^p = G opt , (46) 

^ + c^S - S c n) + S (y + J n N^J +Sx W = ^ G opt o z , (4c) 

— ^ + c n (S c n - Ntf) + — - c' n N S + c'pS c + S c x ft = , {Ad) 
dt t sc 

^ + c' n N n + c' p N lP = 0. (4e) 

Here Pj = (G^ pt — G° pt )/G fopt , o z is the unit vector directed along the normal z coinciding 
with the exciting-light propagation direction, uj and O are the Larmor frequencies 
defined by hw = gfisB, hQ = g c /j,BB, g and g c are the g-factors of free and bound 
electrons, B is the external magnetic field and /is is the Bohr magneton. In the following 
we assume B || y _L z. Equations (14"q|) - (14~e|) together with equations ([T]) and ([3]) form a 
complete set in the proposed model. 

2.2. Continuous wave photoexcitation 

Under continuous wave photoexcitation the time derivatives vanish and the equations 
for S and S c are reduced to 

f -^ + 5x W = ^G opt o 2 , (5) 
T c T 

where the four introduced times are defined by 
1111 

™ = — + - , — = c n n, (6) 

T T s t r c 

7f = \-c n n + c' p p, - = c n N 1 + c' n N . 

By using equations (|S]) we can present the free- and bound-electron spin components 
5* 2 , S cz in the form 

p.Qoptrp 

S *=^M~> S CZ = AS Z . (7) 



where 



9 ~ 9 1 - coVlTT c ) 

* = 1 + cu 2 T 2 - ri- ^- , 8 

1 + fi 2 T 2 

T r 1 — uQTT r 

A = — ^ , 

r l + f2 2 T 2 
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and 7] = (TT c /tt c ), T = T/ \/l — rj, T c = T c j y'l — r\. The transverse spin components 
S x , S cx (for B || y) are given by 

T ujt c S z + QT C S CZ T c uTS z + Vtr S cz 

T c 1—7] T 1 — T] 

The densities p, N , N 2 can be expressed via n and Ni as 
p =<^ N c -N 1+ n 

4^ 1 N c - iVi + n 



c^n" 21 + «/24)(n/iV 1 ) ' 
= KiVc-iVOll + K/g^/iV^-n 
2 2 l + C^^Cn/JVk) ' 1 J 

Excluding p,N ,N 2 , S and <S C from the complete set of equations we obtain, instead 
of equations (14 a\) and fj46|) . the following two coupled nonlinear equations for n and Ni 

ir (iv.-ivoti + K/c;)^/^)]-^ 



where 



(13) 



x 7 ^T/Vr = G° pt , (12) 
(24/0(^/^ + 1 1 J 

c^A^l - E) + c' n N {n, Njn = G opt , 

_ ASS C _ (P^T) 2 
nN x nN^ 2 

x (a + - V ,J u 2 t c T + AuQtt c (1 +rj) + A 2 Q 2 tT, 

According to equation (J7|) the degree of free-electron spin polarization is given by 

p= 2_S z = P 1 G^T 
n 

The intensity / and degree of circular polarization of the photoluminescence p can be 
found from 

/+ - I 

I = I + + I~ ocnp, p = = P'P , (15) 

where 1^ are the intensities of the o~ ± circularly-polarized components of the interband 
PL, and P' = P[P 2 , with the factor P[ taking into account the selection rules for the 
room-temperature recombination of the conduction electrons with heavy and light holes, 
and P 2 < 1 being the depolarizing factor arising due to possible multiple reflections from 
the sample boundaries [35] . 



2.3. Two- charge- state model 

As follows from references [TBI [T5J [HI [T71 [TBI EH EE] and demonstrated below the SDR 
model taking into account only the deep-center states D 1 and D 2 satisfactorily describes 
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the experiments on spin dynamics in GaAsN samples. Therefore, in what follows we 
focus on this particular model setting d = in which case N Q = 0, N c = N± + N 2 , 
p = n + N c — Ni and equations ( 1T21) , ( ITBl) reduce to 

c p (N c - Ni)(N c -N 1 + n)= G op \ c^nN^l - E) = G opt . (16) 

Optical orientation of electron spins at zero magnetic field. Let us consider 
particular limiting cases allowing for simple analytical expressions. In the absence of 
zero magnetic field the Larmor frequencies oo, Q vanish and the expressions for S z and 
E reduce to 

p. g o P t r (P t G°pt T ) 2 T c 1 

* 2(1 - v ) ' nN x r (1 - V ) 2 ' 1 ] 

It is convenient to introduce the dimensionless variables 

x = ££, y=»^*=%, z = £. (is, 

In terms of these variables one has 

p = N c (Y + Z), IcxZ(Y + Z), p = ff^I-^, (19) 



1-1^ Z 



T 

A = — , * = 1 -77, 
r 

1 1-Y 1 1 1-Y 



r r* T r s r* 

1 Z 1 1 Z 

r c r* ' T c r sc r* ' 

where r* = (c n A r c )" 1 and = (cpAc)" 1 . 

According to equations (fTBl and (fTBT) the variables Y and Z satisfy the equations 



Y(Y + Z)=X, (20) 

F U - P 2 I^Ya 2 g>lV!££ . 

T*J [Z + t*/t sc + (1-Y)t c /t sc ]? ' 

where a = c p /c n . 

In the low-power limit we have r c — >■ 00, 77 — >■ and 

n = r*G opt , p = P'Pi Ts , (21) 



pocN c -N l0 c y/G^ , I oc (G opt ) 3 / 2 . 

The first nonvanishing corrections to n and p are taken into account in equation (|2ip 
by the replacement r* — > r*(l + a/ r^G opt /N c ). In particular, this means that 



However, as soon as a value of (t s t sc /t*t^)\/ X exceeds unity the main correction to p 
originates from the spin-dependent recombination, it is linear in G opt and given by 



t s + t* V T s + T* N c J 



(23) 
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With increasing pump intensity, Y monotonously increases and saturates up to some 
limiting value controlled by the ratio c p /c n . This limiting value = Y(G opt —> oo) 
satisfies a third-order algebraic equation 

(l-Y oo )(l-bY*)-aY oo = 7 (24) 

where b = Pi(T s /rl) 2 . Note that according to equation ([20]) one has asymptotically 
Z(G opt -> oo) ->• X/Foo. If a = c p /c n < |1 - Vb\, then 

f 1 - t^t for b < 1 , 

.(i-t+^fcr6>i. (25) 



%/b V 2(v / 6-l) 

For b « 1 and aCl, one has = 1 — y/a/2. 

Hanle effect at transverse magnetic field. Assuming T ^T c and comparable values 
of g and g c we can separately analyze two regions of the transverse magnetic field: (I) the 
low-field region where uT <C 1 and the product QT C is arbitrary, and (II) the high-field 
region where QT ^> 1 and the product uT is arbitrary. 

In the region I, at low pumping level, equation (123|) derived at B = transfers to 

*B) = P'P,^ (l + _|_^_^) . (26) 

One can see that in this particular case p(-B) is a sum of a constant and a weak narrow 
Lorentzian. 

In the region II, at arbitrary pumping level, one can neglect the influence of bound 
electrons upon the formation of the free-electron density and polarization, set S to zero 
and use the equations 

n = r^, ,(B) = f vp,_^__i_ = i + ( ^ )a , (27) 

where po = P'PiTsjijs + t) and £?i/2 = h/(gfiBT) is the half-width at half-maximum of 
the Lorentzian. At weak incident- light intensities, the lifetime r can be replace by the 
characteristic time r*. 

Electron density for unpolarized carriers. The system is spin-unpolarized under 
linearly-polarized photoexcitation as well as for a pumping in the limit of very strong 
transverse magnetic fields. In the both cases, in order to find n and N\ one can formally 
set Pi = and reduce a pair of equations (120|) to 

Y 2 (l - Y) 

— -4= = X . (28) 

l-(l + a)F V ; 

To trace the saturation of Y and asymptotic behavior of Z with the increasing pump 
intensity, we set X — > oo in equation (128j) and obtain 

N, 



(1 + a) 3 r*G°pt 



(29) 



Z(G^ -> oo) 



X c (1 + a) 
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3. Selection of model parameter values 

The two-charge-state model is characterized by two free (controllable) variables, the total 
generation rate G opt and the strength of magnetic field B, and nine fitting parameters, 
namely, Lande factors g and g c , the spin relaxation times r s and t sc , the recombination 
times r* = (cnNc)' 1 and = (c p iV c ) _1 , the degrees of polarization Pi, P' and the density 
of deep centers N c . In the present section we will estimate values of these parameters 
by using available experimental data and reasonable theoretical considerations. 

1. Lande g-f actors of free and bound electrons 

The recent studies [211 121] of ODMR have provided unambiguous experimental evidence 
that, in dilute nitrides GaAsi-^Na, and Gai-ylnyAsi-^N^ with a few percent of nitrogen 
(x < 3.3%), the deep paramagnetic centers are formed by Ga^ self-interstitials. In 
the above notations, the defect state D\ occupied by a single electron represents the 
doubly positive charged level Ga^ + while the singlet D 2 occupied by a pair of electrons is 
identified as the singly positive state G&] + . The bound-electron g factor g c determined 
from the ODMR measurements has been found to be very close to 2. The positive sign 
of g c has been independently determined from asymmetry in the depolarization of edge 
photoluminescence in a transverse magnetic field (Hanle effect) at oblique incidence of 
the exciting radiation and oblique-angle detection of the luminescence [131 EH]- The 
results of references [131 EU EH |37] are in agreement with the existing expectation 
that the g factors of electrons bound to deep levels, as a rule, exhibit only a small 
renormalization with respect to the gyromagnetic factor g = 2.0023 of a free electron 
in vacuum. 

On the other hand, the g values of free carriers, electrons and holes, in 
semiconductors with small and moderate band gaps can drastically differ from go. 
Particularly, for conduction-band electrons in bulk GaAs the g factor equals —0.44 
at helium temperature. In GaAsi-^Nz with a few percent of N the value of electron 
g factor shifts upwards. The experimental study [381 EE] of the Zeeman splitting 
of <t + and o~~ circularly polarized components of the photoluminescence spectra of 
GaAsi-^Nz (x < 0.6%) have shown that the conduction-electron g factor exhibits a 
sign reversal from negative to positive at x ~ 0.04% and increases abruptly up to 
0.7 in a very narrow compositional window between x = 0.04% and x = 0.1%. For 
0.1% < x < 0.6%, the dependence g(x) has a not well-defined behaviour and fluctuates 
around 0.7 for the highest x values. The striking difference from the value —0.44 is 
attributed to the resonant defect level in the conduction band introduced by nitrogen and 
the repulsion between the GaAs host matrix conduction band edge and a higher lying 
band of localized N resonant states (the so-called band-anticrossing model introduced by 
Shan, Walukiewicz et al [40], see also [H]). The agreement between experimental data 
and theoretical calculation is remarkably improved in the generalized band-anticrossing 
model which takes into account the interaction of the conduction-band minimum with 
localized resonant states formed by nitrogen-atom pairs and clusters [38], [39] . 
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The positive sign of g in GaAsi_ x N x film with x = 2.1% was found at 
room temperature from the Hanle effect asymmetry under oblique excitation and 
detection [37] while its absolute value \g\ rs 0.9 follows from the electron spin quantum 
beats in photoluminescence recorded in a transverse magnetic field after the pulse 
excitation [IB] . The data are in agreement with the value g ~ +1 estimated in 
reference [37] from the relationship between g and electron effective mass m* « 0.14mo 
measured for x = 2% in reference [12]. The positive free-electron g-factor g = 0.97 
was measured in GaAsi-^N^/GaAs quantum well with x = 1.5% by time resolved Kerr 
rotation at 220 K [13]. On the other hand, in reference [H] only a slow decrease of \g\ 
from 0.28 to 0.22 for x increasing from 0.09% to 0.9% has recently been reported from 
time-resolved Kerr rotation measurements in GaAsi-^N^. films at room temperature. 
The difference with the above values may be related to variations in nitrogen-atom 
spatial distribution inside dilute nitride samples studied by the two groups. 

Below, we will use g = +1, g c = +2. 

2. Characteristic free electron and hole lifetimes r* and t)* 

According to equation ( l46|) . in the two-charge-state model (c' p = 0) the hole 
recombination rate equals c p N 2 p. In the limit of high pump intensities, G opt — > 00, 
the hole lifetime Th = (cpA^) -1 tends to a limit of = (c p N c )~ l . The regime of 
high pumping was realized by Kalevich et al [16] under pulsed optical excitation of 
GaAso.979No.021 alloy at room temperature. The ratio of the photoelectron density n 
and N c exceeded 10 and the photoluminescence decay time at the initial stage, (15.0 
± 0.5) ps, could be identified with r^/2. Thus, in this particular sample, the time 
can be taken equal to 30 ps or a little bit shorter because, as the analysis shows, in the 
experimental conditions a maximum value of N 2 was slightly smaller than N c . 

In the low intensity limit, where r — > r*, one can use equations ( J2T1) . Moreover, it 
follows from the experiment [16] that the three characteristic times r*, and t s satisfy 
the hierarchy: r* <C T* h <C t s . Therefore, at low pumping levels the lifetime, T, of 
free-electron spin polarization given by the first equation ([6]) is close to r*. This allowed 
us to determine t* from measurements of the Hanle effect at low incident light intensity 
J. The interpolation of the experimental dependence Bi/ 2 (J) ~ h/(gfiBr) of the free- 
electron Hanle curve half-width to the limit J — > yields [36], for g — 1, the conduction 
electron lifetime r* = t(J —¥ 0) = 1.9 ps. 

For the model calculations we will use r* = 1.9 ps and = 27 ps. 

3. Spin relaxation times of free and bound electrons 

The spin relaxation time r s of free photoelectrons was found [151 EH] from the decay 
curves of 

SM x ' + v -/-w K ^w-rw 

P(*) yW) 
and, independently, from the PL intensity decay, I{t), at the second stage of the 
spin-dependent electron dynamics, the partial intensities 1^ being introduced in 
equation (floT) . The obtained two values (150 ± 15) ps and (144 ± 4) ps, respectively, 
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coincide within the uncertainty limits, confirming the interpretation of the experimental 
data. The spin-Kerr rotation measurements jH] performed on three GaAsN samples at 
room temperature give for the free-electron spin lifetime T the values of 60, 120 and 
125 ps which can serve as a lower bound for the spin relaxation time t s in these samples 
in agreement with the above two values of r s . 

The spin-relaxation time of bound electrons, t sc , can be deduced from measurements 
of the Hanle depolarization curve [131 EE]) or ) to be more precise, the narrow part 
of this curve related to the depolarization of bound electrons. Indeed, according to 
equation (1261) the half- width of the narrow Lorentzian is given by By 2 = h~/g c pBT c 
which reduces, at J — > 0, to h/g c pBT~ sc - By using the values B^ 2 ~ 85 G and g c = +2, 
we obtained [3E] t sc ~ 700 ps. 

Later on we will use t s = 140 ps, t sc = 700 ps. 

4- Initial photoelectron spin polarization Pi and the factor P' relating p with P 

The determination of the initial spin polarization Pj from experiments on polarized 
photoluminescence can be complicated by the presence of a small additional contribution 
of unclear origin to the secondary emission. For the circularly-polarized incident 
light, this contribution is also polarized and its circular polarization is insensitive 
to the transverse magnetic field up to very high fields. Taking into account the 
additional contribution the measured degree of circular polarization is presented as 
Pcxp = (1 — j)p + JPadd) where p is given by equation ( fl5l) . p a dd is the polarization of the 
additional radiation and j is its fraction in the total intensity of the secondary radiation. 
Assuming j to be small, we can approximate p exp to 

Pcxp = P + Pres , (30) 

where p res = jp a dd is the effective "residual" polarization. By using equation ( T30|) one 
can find Pi from the following equation 

w p (B = 5kG)-p rcs = 
p cxp (B = 0,high J) - p rcs 
Here p exp (-B = 5kG) is the degree of circular polarization measured at the moderate 
value of the transverse magnetic field, 5=5 kG, high enough in order to neglect the 
spin polarization of bound electrons and apply equation (1271) and, on the other hand, 
low enough to satisfy the condition ujT <C 1; p cxp (B = 0,high J) is the polarization 
degree measured at zero magnetic field and at a high photoexcitation level where the 
spin-filter effect results in almost 100% spin polarization of the free photoelectrons, 
P(B = 0,high It follows then from equation fTT3|) that 

pl = Pex P (5 = 0,high J) -p res 

P(B = 0, high J) 1 ' 

Furthermore, according to equation ( l27j) . at the moderate magnetic field one has 
Pcxp(-B = 5kG) = P'PiT s /(r + T s ) + Pres- In the studied GaAsN alloys, at the moderate 
intensities, the free-electron lifetime is much shorter than their spin relaxation time and 
Pex P (-£> = 5kG) reduces to P'Pi + p res which leads to Pi = [p exp (B = 5kG) — p TCS }/P' 
and finally to equation (|5Tj) . 
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5. Density of deep centers N c 

The density of deep centers responsible for the spin-dependent recombination can be 
found by fitting the measured dependencies of the PL intensity and polarization on 
the incident light intensity with computation of the kinetic equations (j4aj) — (|4ejl . The 
inhomogeneous terms in these equations are proportional to the generation rate G opt 
which can be expressed in terms of the excitation photon energy E p ^, the incident light 
intensity J, the light absorption coefficient a and the beam spot area on the sample 
surface, S*. Taking into account that the thickness of the GaAsN layer L is smaller 
than ct~ l and assuming that each photon creates one carrier pair we have 

a J 

G opt = — — . (33) 
E ph S* v ; 

Let us estimate a typical value of G opt in our experiments performed on 
GaAso.979No.021 at room temperature and the energy E p ^ = 1.312 eV of the exciting 
photons. For this nitrogen content the room-temperature band gap of the alloy equals 
Eg « 1.106 eV [12], and, therefore, the excess photon energy E p h — E g — 206 meV. We 
take a = 2 x 10 4 cm -1 for this difference E p h — E g in accordance to values 2.4 x 10 4 cm -1 
and 1.6 x 10 4 cm _1 measured in the alloys with x = 5% [IS] and x = 2.3% [16] in the 
wide range of photon energies. For the light spot diameter 0.2 mm and the intensity 
J = lmW the generation rate G opt is 3 x 10 23 cm _3 s _1 . 

A lower estimate for the critical value of the generation rate G^ circ at which the 
spin filter effect becomes remarkable can be obtained by setting E(G opt ) to unity. At 
zero magnetic field, small intensities, T c ps t sc and N\ ~ N c one has from equations ( T131) . 

m 

v _ (R t G opt T) 2 T c _ p2 G opt T 2 T c ^ p2 G opt r sc 



nN\ r(l - t]) 2 1 N c t* 2 1 N c 
and, thus, the critical value of G opt is given by 

N 

^cr^irc p2 ' V ^/ 
i 'sc 

It follows from experimental data presented in figure [2b that the spin filter effect develops 
to full scale at J = 75 mW, so that G° r pt circ ~ 3 x lO^cm-V 1 . Taking P { = 0.24, 
t sc = 700 ps and T = r* (since t s ^> r*) we obtain from equation (1341) the estimation 
N c = 10 15 cm" 3 . The best agreement with experiment, see the details below, is obtained 
for the density of deep paramagnetic centers JV C = 3 x 10 15 cm -3 , in agreement with 
the crude estimation. It is worth to mention that, according to equation (f28l) . the 
critical generation rate G° p \ iD at which the system exhibits nonlinear behaviour under 
linearly polarized photoexcitation is given by the condition X = 1 with X defined by 
equation ( 1181) or, equivalently, by 

= ^ • (35) 
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50 100 150 200 

J (mW) 

Figure 2. Spin-dependent recombination ratios (a) and the photoluminescence 
circular polarization (b) in GaAso.979No.021 vs. the excitation power for the incident- 
photon energy hu cxc = 1.312 eV, the detection energy fiwdct = 1.17 eV and room 
temperature. Full and open squares stand for the measured values of -Ksdr = 
7(circ)//(lin) and -^Qdr 5 = I(circ,B = 0)/I(civc,B = 5kG), respectively. Circles 
and triangles present values of the degree of circular polarization measured in the 
absence of magnetic field and at B = 5 kG. Solid and dotted lines are the result of the 
model calculation (see text for details). 

4. Comparison with experiment and discussion 

One of the most convincing signatures of spin-dependent recombination is dependence 
of the PL intensity on the light polarization at zero magnetic field, or more precisely, 
an increase in the PL intensity under the switching from linear to circular polarization 
of the exciting light. In figure |2] full squares show the measured SDR ratio defined 
by K SDR = I(circ,B = 0)/I(lin,£ = 0), where I{citc,B = 0) and I{lin,B = 0) 
are the PL intensity under circularly- and linearly-polarized photoexcitation. At very 
weak pumping power, the intensities J(circ) and /(lin) coincide. With increasing the 
power the PL intensity /(circ, B — 0) is enhanced as compared with /(lin, B = 0) 
and, at J > 100 mW, the SDR ratio reaches values as high as 8 and even higher. 
A transverse magnetic field of 5kG strength eliminates the enhancement caused by 
the circular polarization of the incident light and, as a result, the modified SDR ratio 
-^SDR 5kG ' > = ^(circ, B = 0)//(circ,i? = 5kG) behaves as a function of J identically 
to the dependence Ksbr(J)- Simultaneously the degree of PL circular polarization 
induced under circularly-polarized photoexcitation (full circles in figure |2]) monotonously 
increases up to 25% which is another important manifestation of the spin-dependent 
recombination. 
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Figure 3. The photoluminescence intensity (a, c) and circular polarization (b, c) 
as functions of a transverse magnetic field measured in GaAso.979No.021 for different 
powers of the circularly-polarized pumping. Experimental points presented in panel 
(a) are obtained for the following incident-light intensity J (in mW): A - 25, V - 
50, o - 100, * - 150, • - 220. Panel (b) shows the upper part of the Hanle curve at 
J = 220 mW and panel (c) shows the dependencies 1(B) (solid circles) and p{B) (open 
circles) measured at J — 220 mW in the wide range of magnetic fields. Solid line in 
panel (c) is the Lorentzian p(B) = po/(l + B 2 /B^, 2 ) + p res with B\/ 2 = 25000 G. 
Swexc = 1.312 cV, tvjJdct = 1.17eV. 



Figure [3] shows typical curves of magnetic-field-induced quenching of the 
photoluminescence (panels (a) and (c)) and its depolarization (panels (b) and (c)) in 
GaAso.979No.021 measured under circularly-polarized pumping. One can see that the 
depolarization curve (Hanle effect) is composed of narrow and wide contours with the 
half-widths differing by two orders of magnitude. The main variation of narrow Hanle 
contour takes place at the same magnetic field as for the switching off the PL intensity 
enhancement which is a clear evidence that the narrow and wide parts of the Hanle 
curve are related to the spin depolarization of bound and free electrons, respectively. 

For a quantitative description of the data presented in figures |2] and |3] we apply the 
two-charge-state model for the selected set of parameter values, i.e., 

g = l, g c = 2, r* = {cnN,)- 1 = 1.9ps, (36) 
t* = (c p N c y l = 27 ps, r s = 140 ps, r sc = 700 ps . 

Three remaining parameters of the theory, Pj, P' and N c , are determined from 
comparison with the experimental data displayed in figures [2] and [3j According to 
equations (J30l and (J32l) values of and P' can be found as soon as the effective 
"residual" polarization p res is estimated. Unfortunately, the range of used magnetic 
field is insufficient to determine p rcs unambiguously: its value clearly lies in the 
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Figure 4. Calculated Hanle curve for the free-electron spin polarization p in a coupled 
spin system of free and bound photoelectrons. Inset shows the larger magnetic-field 
range. Solid curve represents the exact calculation while dashed curve is obtained 
assuming unpolarized bound electrons in which case the spin filter is switched off. 



interval between and 3.3%. Taking into account that, as follows from figure EJ^b), 
Pex P (-B = 0,high J) « 25% and p cxp (B = 5kG) ~ 6.7% and using equations (13~0]) . 
032]) we obtain P' « 0.25 and P { « 0.27 (if p rcs = 0) and P' « 0.22 and P « 0.16 
(if p res = 3.3%). An additional variation of these values follows from the difference 
of P(B = 0,high J) from unity since the ratio t*/t s is small but finite. The self- 
consistent calculation, see figure El shows that one should substitute into equation (132|) 
the value P(B = 0, high J) = 0.9 so that the values P' = 0.25, P = 0.27 are replaced 
by P' = 0.28, Pi = 0.24. Due to the same reason, the values of P', Pi corresponding to 
p res = 3.3% should also be changed from 0.22,0.16 to 0.27,0.14. The curves in figure [2] 
are calculated for the two pairs of P', Pi values with N c being the final fitting parameter. 
One can see that the set 

Pi = 0.24, P' = 0.28 and iV c = 3 x 10 15 cm~ 3 (37) 

rather well simulate the experimental dependencies. 

The complete set ( )36|) . ( 137]) of parameters allows us to simulate the effect of 
transverse magnetic field and analyze consequences of variation of each parameter. The 
theoretical Hanle curve calculated at pump power J = 90 mW is shown in figure H] in 
the narrow range up to 700 G and, in inset, in the wide range up to 60 kG. As mentioned 
in section El the Hanle curve outside the narrow contour is described by the Lorentzian 
( |27|) . This property is demonstrated by drawing a subsidiary (dashed) line calculated 
assuming formally that paramagnetic centers are kept unpolarized, i.e., setting S c = 
in equations ( I4~ffl|) - (j4~e|) . Then the densities n, N\ and p are independent of the magnetic 
field and the polarization p(B) is indeed given by the Lorentzian ( |27l) . In figure H] the 
exact (solid) curve is rather well approximated by the dashed curve for B > 0.5 kG. 

The two-charge-state model rather well reproduces manifestation of the Hanle effect 
in the PL intensity. The dependence of PL intensity on the transverse magnetic field 
calculated for different a-pump intensities is shown in figure [5j In agreement with 
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Figure 5. Magnetic- field dependence of the PL intensity 1(B) calculated at different 
powers J of circularly polarized excitation and presented on small (a) and large (b) 
field scales. The curves 1. . .5 correspond to the following values of J: 1 — 10 mW, 
2 - 50 mW, 3 - 90 mW, 4 - 170 mW, 5 - 220 mW. 



experiment the intensity / decreases in relatively weak fields | B \ < 5 kG which suppress 
the polarization of bound electrons. Moreover, half-width of the decaying part of 
the function 1(B) is close to that of the narrow contour of the curve p(B). At high 
fields \B\ > 5kG gradually depolarizing the free electrons, the PL intensity is field- 
independent, see figure E^b), in accordance with experimental data of figure Etc). 

Figure [6] displays the relation between the average electron spin polarization P 
and the initial polarization Pi for different pumping powers J. At low values of J this 
relation is described by equation (1231) which can be rewritten in the form 

p= { c " + Cl T^w) P - (38) 

where c = T/r* and Ci = clT sc G opt (lmW)/N c . Thus, at weak pumping and small Pi, 
the spin polarization P is proportional to Pi with the proportionality coefficient being 
a linear function of J. This approximate relation is shown in figure El by dotted straight 
lines. One can see that, at small J, the onset slopes of the exact and approximate curves 
P(Pi) coincide; at moderate pumping the slopes slightly deviate, and only at large values 
of J the difference between slopes gets remarkable. With increasing Pi the dependence 
P(Pi) at finite fixed J becomes sublinear and finally saturates to r s /(r + t s ) m 1. 

Figure [7] illustrates effect of the sign of free-electron g factor on the Hanle curve. The 
curves in the left and right panels are calculated for g = 1 and g = — 1 , respectively, while 
the parameter g c = 2 is kept constant. One can see that, the value of p(B = 5 kG) (which 
determines the amplitude of the wide Hanle contour) calculated at g = 1 remarkably 
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Figure 6. Relation between average spin polarization, P, of the conduction 
photoclcctrons and their initial spin polarization Pi calculated for different excitation 
powers varying from lOmW to 400 mW. The dashed line shows the limiting relation 
P = PiT s /(r + t s ) ps Pi attainable at vanishingly small intensity where the spin filter 
is switched off. The dotted straight lines represent the linear approximation ([38]) . The 
shaded area Pi > 0.5 is unachievable in a bulk undeformed semiconductor. 

decreases with the increasing J. In contrast, at g = — 1 the wide contour moves upwards 
as the intensity J increases from 10 mW to 250 mW and, moreover, the narrow contour 
becomes much wider as compared to the case of g = 1. This can be understood taking 
into account that, according to equation (jHJ), in the system of coupled free and bound 
electron spins they precess in the transverse magnetic field faster if the signs of g and 
g c coincide and slower if the signs of g and g c are opposite. As a result, for negative g 
and positive g c the function \1/ in equation (jSJ) becomes more robust to the influence of 
the transverse magnetic fields in the region B < 2kG. Comparison of the experimental 
Hanle curves in figure 2 of reference [36] with figures 0(b) and (d) confirms the positive 
sign of g determined in GaAso.979No.021 under oblique pumping and detection [37] . 

The behaviour of the system under study in the wide range of pump intensities 
is demonstrated in figure [HJ At weak intensities the dependence of N\ and N 2 is 
insensitive to the polarization of incident light: Y = A^/Ac = (N c — N\)/N c oc \/~J. 
With increasing J the dependence A 2 (J) tends to saturation. For r* = 1.9 ps, t sc 
= 700 ps and Pj = 0.24, we have for the critical values ( )34|) . (|35|) of the generation 
rates: lin ^> G°^ circ . Therefore, the deviation from the square- root law begins earlier 
for the circular polarization. The limiting values of Y at G opt — > 00 are the solution 
Yoo(circ) = 0.72 of equations ( |24l) and the solution Y^flm) = c n /(c n + c p ) = 0.93 of 
equation ( 128]) at X — > 00. The limiting values of Asdr an d P are [(l+a)F 00 (circ)] -2 = 1.7 
and (P / P i r sc /r^)y^ (circ) = 0.26, respectively. The free-electron spin polarization tends 
to Poo = (pr sc /r^)Y 00 (circ) = 0.9. Note, however, that these limiting values are 
achieved at very high intensities far exceeding J = 1000 mW where the model defined 
by the parameter set ( |36|) . ( 1371) is scarcely applicable. Particularly, an increase in 
the photoelectron and photohole densities at very high generation rates can lead to 
a remarkable contribution of the interband electron-hole recombination and the term 
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Figure 7. Hanle curves calculated for positive, g = +1 (a, b), and negative, g = — 1 
(c, d), values of the conduction electron g factor and for different excitation powers 
J = 10 mW (dashed), 90 mW (dotted) and 250 mW (solid). Other parameters are 
indicated in text, see (l36l) and (1371). 



^n±p in equation (EJ) should be taken into consideration. A nonmonotonic behavior 
of -Ksdr as a function of J is an interesting and intuitively unexpected. In figure [2] 
experimental points only begin to show such behaviour. However, references [20l E3] 
report observation of a maximum in the dependence Ksbr{J)- As one can see from 
figure El(b), the degree of PL circular polarization rapidly increases with increasing J 
and reaches a maximum value p max = 0.27, then it almost stops to depend on J slowly 
decreasing to = 0.26. The corresponding values of P are P max = 0.94 and Poo= 0.90. 



5. Conclusion 

In summary, we have developed a spin- dependent theory of Shockley-Read 
recombination for three- and two-charge-state defects. The two-charge-state model 
characterized by nine parameters satisfactorily describes the coupled spin system of free 
and bound photoelectrons in undoped GaAsN alloys. In the weak- and strong-pumping 
limits, we present simple analytic equations which allow us to predict many properties 
of the coupled system not applying to computer simulations. Here we have concentrated 
on cw optical excitation. The further step could be a comprehensive analysis of time- 
resolved kinetics, including PL decay after short-pulse photoexcitation and pump-probe 
spectroscopy. It might be also interesting to study the strong-pumping regime where 
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Figure 8. Variation of Kson, p, Ni/N c , N 2 /N c , N+/N c and N-/N c with the 
increasing pump intensity J calculated in the wide range of J for the set of parameters 
(f36f and (f37]) . Short horizontal lines indicate the limiting values at J — >• oo. 



the recombination of free electrons and holes acquires an important role in the kinetics 
of photocarriers and work out predictions of specific behaviour of the system in this 
regime. Experimentally, it would be rather instructive to demonstrate unambiguously 
an existence or nonexistence of the deep-center level in GaAsN with no electron on the 
center. 
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